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Abstract 

. . . A long standing problem in Lattice QCD has been the discrepancy between the experimental and calculated values for the axial 
charge of the nucleon, gA = Ga(Q^ - 0). Though finite volume effects have been shown to be large, it has also been suggested that 
excited state effects may also play a significant role in suppressing the value of gA- In this work, we apply a variational method to 
generate operators that couple predominantly to the ground state, thus systematically removing excited state contamination from 
^ the extraction of gA- The utility and success of this approach is manifest in the early onset of a clear plateau in the correlation 
function ratio proportional to gA- Through a comparison with results obtained via traditional methods, we show how excited state 
effects can suppress gA by as much as 12% if sources are not properly tuned. 
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> 1. Introduction 

In recent years, lattice calculations have taken a tremendous 
step towards simulating QCD at the physical point. Algorith- 
mic and technological developments have allowed simulations 
to probe at or near physical quark masses on increasingly larger 
Volumes, with finer lattice spacings and vastly increased statis- 
tics. Calculations of the ground state spectrum have yielded 
results consistent to within a few percent of their physical val- 
ues with well controlled systematic errors [1, 2]. Naturally the 
next step has been to strive for this level of precision for the 
jnatrix elements of these states. Despite the remarkable consis- 
tency between lattice and experimental data for the pion form 
factor F„(Q^), a complete description of other hadronic states, 
particularly the nucleon, has proven to be remarkably challeng- 
^ jng [3,4]. 

. ^ ' The most notable shortfall is for the nucleon axial charge, 
gA = Ga(Q^ - 0). In principle should be relatively simple to 
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Calculate. Being an iso-vector quantity, disconnected loop con- 
•tributions are absent and as we have direct access to Ga{0), we 
circumvent the need for extrapolations in Q^. Unfortunately, 
the lattice values for gA to date have been consistently lower 
than the experimental value by as much as 10-15% [5]. In an 
effort to account for these discrepancies, several studies have 
carefully examined the systematic errors present in the calcula- 
tion [6-16]. In this letter we will focus on the role of excited 
state effects. 

Recently there has been an increased effort to understand and 
reduce the impact of excited states on form factor calculations. 
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In computing these quantities, it is well understood that to en- 
sure excited state contributions to the correlation function are 
sufficiently suppressed, one needs large Euclidean time separa- 
tions between operators. To choose a suitable time separation 
one should identify the time slices where the correlation func- 
tions take on their asymptotic form. For the two commonly 
used sequential source techniques, this is a relatively simple 
procedure for the fixed current method. One simply chooses a 
current insertion time, tc, once the asymptotic behaviour is ob- 
served in the two-point correlator Results are extracted from 
the data once the asymptotic behaviour is observed in the three- 
point correlator. 

For the fixed sink method, one requires knowledge of the 
asymptotic behaviour of the three-point correlator a priori- Un- 
fortunately, the temptation to use earlier sink times in order to 
obtain more precise results is inescapable. These results can 
suffer from excited state contaminations, even if a plateau is 
observed with tc- In ref. [14], it was found that for certain ma- 
trix elements, eg. (x), the source-sink separations often used 
in the literature were not sufficiently large to suppress excited 
state effects.. Nonetheless, as we move ever closer to the physi- 
cal point we are naturally forced to choose earlier sink times as 
signal degrades much quicker 

To counter this issue, new techniques are being devised to 
try and control the sub-leading terms to the three point corre- 
lator The use of the summation method [16, 17] has shown 
improvement upon the conventional approach, but the underly- 
ing excited states contributions are still present. It is not hard 
to imagine situations where these still impact significantly and 
alter the final result. 

In this paper we take a somewhat different approach. Rather 
than reduce the impact of excited states through Euclidean time 
evolution, we seek to separate them out from the ground state 
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at the source and sink. Drawing upon the techniques devel- 
oped for excited state spectroscopy calculations, we will use the 
variational approach to construct interpolating fields that couple 
with individual energy eigenstates and use these to isolate the 
desired matrix elements [18, 19]. An analogous approach has 
been presented in [20, 21] for the study of B* — > Bn transitions 
and proved to be remarkably successful. Here we apply it to gA- 
This letter is organized as follows. In Section 2 we will 
examine the variational method in the context of excited state 
spectroscopy and then outline how this method can be applied 
to the calculation of hadronic matrix elements. Section 3 out- 
lines the details of this calculation. In Section 4 we present our 
results and compare our variational method with the traditional, 
single-operator approach to the calculation of gA- Finally we 
provide our concluding remarks in Section 5. 

2. Variational Method for Matrix Elements 

The 'variational method' [22, 23] is a well established ap- 
proach for determining the excited state hadron spectrum. It 
is based on the creation of a matrix of correlation functions in 
which diff'erent superpositions of excited state contributions are 
Unearly combined to isolate the energy eigenstates. A diversity 
of excited state superpositions is central to the success of this 
method. 

Starting from a basis of operators {Xiix) \i= I,.. .,N}, we 
construct a correlation matrix of two-point correlation func- 
tions. Due to the discrete nature of the lattice, we can decom- 
pose these correlation functions into a discrete sum over energy 
eigenstates, 

G,iP,n r) = 2.-^^'zr(p-)Zj(p-)tr (^1^) , (1) 

where the parameters Z?(^ are the coupling strengths of the 

interpolators with the energy eigenstate of mass nia and F 
projects out the desired parity. We choose new operators (p^ix) 
to be Unear combinations 



(2) 



with a suitable choice of coefficients vf and u", such that these 
interpolators couple to a single energy eigenstate. 



{n\<(/'iO)\a,p,s} = d'"'Zaip) 



u(p, s) . (3) 



From Eqs. (1) and (3) we find that the necessary values for v? 
and u"^ are the solutions of the following eigenvalue equations 

vm to + AO {G{p, to)r% = c" v%p) , (4) 

[{Gip, ro))-i G{p, to + mij u"j{p = c" ufip , (5) 

where the eigenvalue c" = e'""^*. 

It is important to note that both (4) and (5) are evaluated for 
a given momentum p and so the diagonaUsation condition is 



only satisfied when we project with the relevant coefficients as 
foUows: 

v''i{p)Gi0,t;T)tf.{p)^6"P . (6) 
Thus the two-point correlation function for the state | a, p) is 

G^p, t; F) = v'^ip Gijip t; F) u'](p . (7) 

We can extract ffie mass ot„ from G"(p = 0, t) in ffie standard 
way. 

To understand how we can utilise ffie variational method for 
use in form factor calculations, we must firstly identify the 
terms present in ffie three-point correlation function. 



Gl'j{p',p t2,h; F') = Yj 

tr ( r (Q \xiix2) J^ixi)xj{0) (8) 

Sandwiching ffie current operator between two complete 

sets of states we end up with three terms, the vertex 
amphtude, {fi,p',s'\J'^(Q)\a,p,s), and ffie coupling terms 
<Q|Ar,(0) I A s'} and (a, p, s \xj{0) \ Q), 



G^jiP',p; h,h; F') = ^ g-^^^''^-''' e"^"^'' 

a,l3 

{l5,P,s'\nO)\a,p,s)u{p,s)y (9) 

The coupling parameters take the same form as they did in the 
calculation of ffie two-point correlator wiffi two key diff'erences. 
The inclusion of a current means that the initial and final mo- 
menta need not be the same. Furthermore, there also exists the 
possibility that the initial and final energy eigenstates are not 
ffie same. That is, the current can induce a transition between 
states. For this calculation ffie necessary expression is 

G'^^C^',^; t2,h; F') = v'{{p)G^.{p',p; HM\ r^u"^{p) . (10) 

To isolate the matrix element from the three-point function, we 
construct a ratio in ffie standard way. In this work we shall use 
the ratio defined in [24]. For ffie state a the necessary ratio is. 



R"^{p',p; F',F) = 



G'^t'ip'J; t2,h; T')G"''{p,p'; tiJu FQ 
' G''{p,t2;T)G»{p',t2;T) 



(11) 

Key to this approach is the utilisation of a basis of operators 
in which there is diversity in the overlap wiffi various excited 
states. As there are a hmited number of local bilinear operators 
for a given Z^*-, a great deal of work has been made by various 
groups in increasing the number of available operators. Here we 
choose to use fermion source and sink smearing as a meffiod of 
extending our operator basis, as outUned in [25, 26]. 
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3. Calculation Details 

For this calculation we make use of the PACS-CS (2+1)- 
flavour dynamical-QCD gauge field configurations [27] made 
available through the ILDG [28]. These configurations are 
generated using a non-perturbatively <9(fl)-improved Wilson 
fermion action and Iwasaki gauge action. The value /3 - 1 .90 
results in a lattice spacing a - 0.091 fm, determined via the 
static quark potential. With dimensions 32^ x 64, these ensem- 
bles correspond to a spatial length of L = 2.9 fm. As the inten- 
tion of this paper is to examine whether the variational approach 
is an improvement upon traditional techniques, we will consider 
only the light quark mass that corresponds to m„ * 290 MeV. 
The resulting value of m„L = 4.26 is comparable to the values 
used by most groups. 

In this work we are primarily interested in isolating the 
ground state and so have chosen to use a small variational basis 
upon which to perform our correlation matrix analysis. We use 
gauge-invariant Gaussian smearing in the spatial dimensions 
only with smearing fraction a - 0.7 [29]. We consider four 
levels of smearing with the optimal choice found in [29], these 
being 16, 35, 100 and 200, applied to the standard, local proton 
interpolator 

thus allowing for construction of a correlation matrix of dimen- 
sion up to 4 X 4. In Table 1 we list the rms-radii for our choice 
of smearing parameters. We choose to use variational parame- 
ters = 18 and = 2, again taken from [29], where it was 
found that this choice produced best balance between system- 
atic and statistical uncertainties. 

Table 1 : The rms radii for the various levels of smearing con- 
sidered in this work. 



Sweeps of smearing 


rms radius (fm) 


16 


0.216 


35 


0.319 


100 


0.539 


200 


0.778 



To extract the nucleon axial charge we are interested in the 
matrix element ( p{p' , s') \ A'lf \ n(p, s) ) where A'^'^ - uy^y^d. 
This vertex can be expressed via two independent form factors, 
the axial form factor Ga(2^) and the induced pseudoscalar form 
factor GpiQ'^), as 



(pip',s')\Af\nip,s)) = y^'j 
Upip', s') 



1/2 



w5GA(e') + r5^Gp(e') 



u„(p,s), (12) 



where q^i - p'^^- Pfi and Q — - q . Using isospin symmetry, 
one can show that the flavour-changing current 

( p(p', s) I A"" I n(p, s)) = { Pip', s) I A"-" I Pip, s) } , 



where A^"'' = uyi^ysu - dyf^ysd. We choose to calculate gA 
using A^""*. 

As we are interested in Ga(6^ = 0), it suffices to consider the 
case where the incoming and outgoing momenta are the same, 
in particular we choose to work in the nucleon rest frame as 
this will provide the smallest statistical uncertainties. This will 
mean that the left and right eigenvectors required to project out 
the three-point function will now correspond to the same mo- 
menta. 

We choose to insert our fermion source at /q = 16. For the 
calculation of the three-point functions we use a local axial vec- 
tor current calculated using a sequential source technique with 
the current held fixed and inserted at = 21, well past the 
onset of asymptotic behaviour for the projected two-point func- 
tion. We choose to use // = 3 for the current with the corre- 
sponding projection matrix being F' = = r4y5y^, where 
F4 = j(7-i-7o)- The value for the axial renormalisation constant 
Za = 0.781 was determined non-perturbatively in [30] using a 
Schrodinger functional scheme. 

The resulting expression from which we extract gA is 



v;'(0)G;^.(0,0;f2,fi;r3)M';(0) 
vO(0)Gy(0,f;F4)M»(0) 



(13) 



As a comparison, we also evaluate gA using a single correla- 
tor from smeared source to point sink and smeared source to 
smeared sink. These are indicative of results one would ex- 
tract from a traditional approach. We examine the four source 
smearing levels used in the correlation matrix calculation. 

4. Results 

In Fig. 1 we present the bare values of gA with increasing sink 
time ts for the smeared source to point sink, smeared source to 
smeared sink (both with 35 sweeps of smearing) and our vari- 
ational method respectively. Between the traditional approach 
(upper two plots) and the variational approach (bottom plot), 
we can see significant difi'erences in the overaU shape of the 
plateau. 

For the smeared source to point sink (upper plot) the Eu- 
clidean time suppression of excited state contributions mani- 
fests itself as a steady increase in the extracted value. This 
trend in the data does not have a clear endpoint and so there 
is no definite plateau. Guided by thex^of obtained via a covari- 
ance matrix analysis, the earliest time slice one could consider 
is ts =25, but what is clear is that we are forced to consider fit 
windows uncomfortably close to regions dominated by noise. 

By smearing the sink as well as the source, there is a definite 
improvement in the quality of the plateau. The excited state 
behaviour is again present as a steady increase in the value of 
gA, but somewhat suppressed. In this case there is a definite 
plateau observed at ts = 24, which is supported by the x\ov 
Unfortunately, this is again somewhat close to the region where 
signal is lost to noise. 

In Fig. 1 (c) we see quite a different situation. Our variational 
approach yields extremely clean results with rapid ground state 
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Figure 1: A comparison of un-renormalized gA as a function 
of sink time. The first two figures are using the traditional ap- 
proach of smeared source point sink and, smeared source 
smeared sink, both for 35 sweeps of smearing. The final figure 
is the result from a 4 x 4 correlation matrix. 
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Figure 2: An overlay of the results from Fig. 1. The data sets 
have been offset from the time slice for clarity - the circles 
(blue) are the results for the variational approach, the trian- 
gles (purple) are the smeared source — > smeared sink, while 
the squares (red) are the smeared source point sink. The fit- 
ted value from the variational approach has been included (blue 
shaded region) to highlight where the traditional approach is 
consistent with the improved method. 

dominance. The systematic rise in the data is no longer present 
and the onset of the plateau is within two time slices of the 
current insertion. 

In Fig. 2 we have overlaid the three datasets to highlight the 
excited state behaviour between the traditional and variational 
approach. If we look carefully at the variational approach, we 
can see that some excited state contamination is present imme- 
diately after the current, but this is short lived. It is worth noting 
that this is a consequence of the limited size of our variational 
basis. As is highlighted in [21], an n x n correlation matrix al- 
lows one to isolate out the n lightest states in the given channel 
and so the sub-leading contributions will come from the n''' + 1 
state. In the case of the ground state, these contributions will be 
short-lived due to the large mass splitting between the ground 
state and n''' + 1 excited state. If one were to construct a basis 
whose dimension was the number of states in the given channel, 
then it would be possible to completely isolate each state. 

What is of most concern in Fig. 2 is the lack of overlap be- 
tween the results of the traditional approach and those of our 
variational method which is free of excited state contamination. 
In Table 2 we list those fits, for the three data sets with the strict 
criterion that the lies between 0.800 and 1.200. In both data 
sets employing the traditional approach, we can obtain good fits 
with small uncertainties if we choose to begin fitting around 
ts - 25 or 26, but find that the results are significantly small. 
The correct result can be extracted from these datasets if we 
choose to fit at later time slices around - 28, but the resulting 
values have unattractively large uncertainties, as they are close 
to the onset of noise. Alarmingly, as we move the fit window 
to later times, the central value increases. Between 25-30 and 
28-30 we observe a systematic variation of 6% in the value g^. 
It is clear that in this case, we have little control over the excited 
state systematics. In contrast, between the various fit windows 
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Figure 3: Comparison of the renormalized value of gA- The first 
four pairs of points are the results for the traditional, point sink 
(squares) and smeared sink (triangles) approach with increasing 
levels of smearing to the right. The rightmost point (circle) is 
the result extracted using variational approach. There is a clear 
dependence on the level of smearing to the extracted result. 



on the data from the variational approach we find the variation 
in the fits is considerably smaller than the smallest statistical 
uncertainty. 

It is worth considering how the level of smearing affects the 
extracted value of gA- In Fig. 3, we present the renormalized 
gA considering each of the four smearings used to construct our 
variational basis. What we find is a dependence on the level of 
smearing used in the calculation. It appears that for low levels 
of smearing the extracted result can be significantly lower, with 
the smallest level of smearing differing by up to 15% from our 
improved, variational result. From this evidence, it is clear that 
if the smearing level is not properly tuned then excited state ef- 
fects significantly impact the extracted result for gA- Increased 
smearing will improve the extracted result, but there is an opti- 
mal value. 

In principle, one could tune the smearing so that the opti- 
mal overlap is observed with the ground state. By using a 
point source propagator and tuning the smearing through the 
sink via the two-point correlator, outlined in [31], one removes 
the need for expensive inversions for each smearing. Unfortu- 
nately, there is no guarantee that the optimal level of smearing 
for a given quark mass, /3 value, momentum or operator wiU 
be universal and so one must tune the smearing for each set of 
parameters under consideration. Immediately, one can find ap- 
peal in the variational approach as there is no longer a need to 
tediously tune the operators to match the ground state. 

The variational approach provides us with a systematic 
framework for constructing operators whereby we have not sup- 
pressed, but instead removed the contributions of the nearby 
states. To see how small one could make the variational basis 
so as to obtain the correct result, we examined all possible sub- 
sets of our variational basis. The results are displayed in Fig. 4. 
To ensure excited state effects are well suppressed it appears 
that the higher levels of smearing are key. Furthermore, clean 
results require at least a 3 x 3 correlation matrix. 
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Figure 4: Results for gA with different number and combina- 
tions of operators used in the variational analysis. 

5. Conclusion 

In this letter we have illustrated how the variational approach 
can be used to eliminate excited state effects from the calcu- 
lation of the nucleon axial-vector coupling constant gA- These 
effects act to suppress lattice simulation results for gA- 

The method is general and would be ideally suited to calcu- 
lations of form factors where the variational approach could be 
applied separately for each choice of source-sink momentum 
combination. Another quantity that has so far proved elusive 
for lattice calculations and could benefit from our approach is 
the quark momentum fraction, (x), which is notorious for pro- 
ducing lattice results that are more than 50% larger than phe- 
nomenological determinations (see [32] for a review). 

Future investigations will accurately calculate gA at a variety 
of quark masses and connect these results to Nature via finite- 
volume chiral effective field theory. 
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